Unidirectional Pinning and Hysteresis of Spatially Discordant Alternans in Cardiac Tissue 
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Spatially discordant alternans is a widely observed pattern of voltage and calcium signals in cardiac tissue 
that can precipitate lethal cardiac arrhythmia. Using spatially coupled iterative maps of the beat-to-beat dynam- 
ics, we explore this pattern's dynamics in the regime of a calcium-dominated period-doubling instability at the 
single cell level. We find a novel nonlinear bifurcation associated with the formation of a discontinuous jump 
in the amplitude of calcium alternans at nodes separating discordant regions. We show that this jump unidi- 
rectionally pins nodes by preventing their motion away from the pacing site following a pacing rate decrease, 
but permitting motion towards this site following a rate increase. This unidirectional pinning leads to strongly 
history-dependent node motion that is strongly arrhythmogenic. 

PACS numbers: 87.19.Hh, ()5.45.-a,89.75.-k 



The study of period-two dynamics in cardiac tissue has 
become an important topic of research in the physics Ol 
and biomedical communities f^. The term alternans de- 
scribes beat-to-beat alternations of both action potential du- 
ration (APD) and peak intracellular calcium concentration 
(CaP^'^'^). Heart cells generically exhibit alternans when they 
are paced rapidly or in pathological conditions. Interest in 
alternans during the last decade has stemmed from the dis- 
covery that APD alternations can become "spatially discor- 
dant" in tissue meaning that APD alternates with oppo- 
site phases in different regions (flEtl- Spatially discordant al- 
ternans (SDA) dynamically creates spatiotemporal dispersion 
of the refractory period during which cells are not excitable, 
thereby promoting wave blocks and the onset of lethal cardiac 
arrhythmias 

To date, our theoretical understanding of SDA is well devel- 
oped for the case where APD alternans results from an insta- 
bility of membrane voltage (Vm) dynamics at the single cell 
level, which originates from the restitution relation between 
APD and the preceding diastolic interval (DI) between two 
action potentials. Numerical simulations [6] have shown that 
"nodes", which are line defects with period- 1 dynamics sep- 
arating discordant regions of period-2 oscillations of opposite 
phases, can form spontaneously in paced homogeneous tissue 
due to conduction velocity (CV) restitution, the relationship 
between action potential propagation speed cv and DI. In ad- 
dition, node formation has been understood theoretically in 
an amplitude equation framework |7|, l3l to result from a pat- 
tern forming linear instability that amplifies spatially periodic 
stationary or traveling modulations of alternans amplitude. 

Despite this progress, our theoretical understanding of SDA 
remains incomplete. Both experiments |9] and ionic model 



SDA formation for Vm- and Cai-driven alternans flliitl^ but 
more complex behaviors for the latter case depending on the 
strength of Cai-driven instability 11211 and the nature of Caj- 
Vm coupling III llllSll . However, a theoretical framework to in- 



terpret both computational and experimental observations has 
remained lacking. 

In this Letter, we extend the theoretical framework of 101 to 
uncover novel aspects of SDA formation for Cai-dominated 
instability and validate our theoretical predictions with de- 
tailed ionic model simulations. A major finding is that node 
motion can be pinned in one direction owing to the formation 
of a discontinuous jump in calcium alternans amplitude at a 
node where only Vm exhibits period- 1 dynamics. This jump 
leads to strongly history-dependent SDA evolution and also 
alters fundamentally the spacing between nodes. We summa- 
rize here our main results and additional details of both theory 
and simulations will be described elsewhere. 

We start our analysis from the system of spatially coupled 
maps of the general form 



Cn + l{x)^ f,[Cr.{x),D„{x)], 

An+l{x)= f Gix,x')fa[C„+l{x'),D„{x')]dx', 
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simulations IIIOII have shown that Caf°^'' can alternate even with G{x) = H^{x) 



where A„{x), D„{x), and Cn{x) denote the APD, DI, and 
^^pcak^ respectively, at beat n and position x, and G{x,x') 
captures the cumulative effect of electrotonic (T4i -diffusive) 
coupling during one beat. For a cable of length L paced at 
X = with no flux boundary conditions on Vm at both ends, 
G{x,x') = G{x ~ x') + G(x + x') + G{2L - a; - x'). 



when Vm is forced to be periodic with a clamped action po- 
tential waveform, demonstrating that alternans can also result 
from an instability of intracellular calcium dynamics. Al- 
though alternans are presently believed to be predominantly 
Cai-driven in many instances, our understanding of nodes in 
this important case remains limited. Numerical simulations 
have shown a qualitatively similar role of CV-restitution in 



1 + ff (l-fr 



where iJ^ is Gaus- 



sian with standard deviation ^ (see Appendix B of 101) and 
£, = \/2-DyAPD* and w = 2Dv / cv* are two intrinsic length- 
scales expressed in terms of the APD and CV at the alter- 
nans bifurcation (APD* and cv*, respectively), and Dy is 
the diffusion constant of Vm in the standard cable equation 
Vm = Dyd^Vm — lion- Furthermore, CV-restitution causes 
the activation interval Tn{x) = A„(.t) + Dn{x) to vary from 
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beat to beat along the cable as QIS 

r dx' 

Tn{x) 



dx' 



cv{Dn{x')) Jq cv{Dn-i{x')) 



(3) 



where r is the imposed period at the paced end (x = 0). 
To complete the model, we need to specify the forms of fa 
and fc. Since we are interested in understanding the generic 
behavior of alternans, we choose simple phenomenological 
forms of those maps defined implicitly by 



fc/ C* = 1 - rc„ + cl + adn, 

fa/ A* = 1 + (idn + 7Cn+l: 



(4) 
(5) 



where c„ = (C„ - C*)/C\ d„ = (Ai - D*)/A*, and we 
also define a„ = (A„ - A*)/ A*. With this choice C„ = C* , 
An = A*, and £)„ = D* are trivial fixed points correspond- 
ing to c„ = a„ = c?„ = 0. Moreover, Cn, a„, and d„ mea- 
sure the departure of Ga^'^^'^, APD, and DI from those fixed 
point values during alternans. The cubic polynomial in c„ in 
Eq. (|4| models a period-doubling bifurcation of the intracellu- 
lar calcium dynamics with the amplitude of Ca?'''''^ alternans 
increasing with the degree of calcium instability r. The term 
dn in Eq. (|5]) incorporates APD-restitution. The other cross 
terms in Eqs. (|4]i and (|5]l model the bi-directional Vm — Cai 
coupling taken to be positive in both directions (a > and 
7 > 0), corresponding to the typical case of locally in-phase 
APD and CaP°^'^ alternans. 

To complete the derivation of maps describing the dynam- 
ics of a„, dn, and c„, we linearize Eq. ((Sj about the fixed 
point, which yields 



dn = - 



J _(d„-d„„i)/2 



(6) 



where A = cu*^/(2cw'*). This linearization is formally valid 
as long as the amplitude of DI alternans induced by Caj alter- 
nans is small enough that we can locally neglect the curvature 
oftheCV-restitutioncurve,or |(r-l)i/27A*c-(;"*/cv'*| < 1. 
Furthermore, we assume that the evolution of alternans ampli- 
tude is sufficiently slow that we can make the approximation 
dn-i ~ —dn- This assumption is valid close to the alternans 
bifurcation. Substituting (i„_i = — (i„inEq. (|6]l and differen- 
tiating both sides, we obtain a differential equation for dn{x) 
that can be solved exactly. Because the pacing rate is fixed at 
X = 0, giving a„(0) + dn{0) = 0, this yields 

dn{x) = -an{x) + e--/^ / /^a„(a;'). (7) 



A 

The dynamics is completely specified by Eq. (|7]) together with 
the maps obtained by inserting Eqs. dUi-llS]) into Eqs. ([T]!-©: 

Cn+i{x) = -rcn{x) + c\[x) + a(i„(x), (8) 

an+i{x)= [ G{x,x')[(3dn{x')+-iCn+i[x')] dx'. (9) 
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FIG. 1 : (Color online) (a) Nature of steady-state solutions in the r-A 
plane. From left to right, no alternans (c = 0), smooth calcium pro- 
files (c > 0, smooth), and discontinuous calcium profiles (c > 0, dis- 
continuous), (b) Smooth traveling and (c) discontinuous stationary 
c{x) (blue) and a(x) (dashed red) profiles from simulating Eqs. (O- 
using (r, A) = (0.9, 10) and (1.2, 15), respectively. Alternans 
profiles obtained from a detailed ionic model jlS ] analogously show- 
ing (d) smooth traveling and (e) discontinuous stationary solutions. 

Note that the pacing rate r no longer appears in the final equa- 
tions but is still contained implicitly in the fact that D* , and 
hence the CV-restitution slope and A, can depend on r. Also, 
since dn-i = — rf„ in steady-state, Eqs. (I7]i-(|9|l remain valid 
in steady-state even further from the bifurcation. 

In Fig.[Tl we present the results of different alternans behav- 
ior obtained from a numerical survey of Eqs. ©-(HI) where we 
vary systematically CV-restitution, which becomes shallower 
with increasing A, and the strength of Cai-driven instability 
that increases with r. In Fig. [TJa) we summarize the nature 
of steady-state solutions in this parameter space. Small r val- 
ues yield no alternans solutions (c = 0) where both steady- 
state c{x) and a{x) are identically zero. When r is increased 
we find a first bifurcation at a value ri (A) where steady-state 
solutions for both c{x) and a{x) become non-zero and form 
smooth waves (c > 0, smooth). If the asymmetry of G, given 
by w, is not too small these waves are stationary, otherwise 
they move towards the pacing site with a constant velocity, 
as in the voltage dominated case ff]. For all work presented 
here w = was used, yielding traveling waves in the smooth 
regime. We found qualitatively similar results for positive w. 

When r is increased further we find a second bifurcation at 
a value r2(A) where calcium alternans profiles become sta- 
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tionary and discontinuous at the nodes separating out of phase 
regions (c > 0, discontinuous) while a{x) remains smooth 
due to the smoothing effect of vohage diffusion. Example 
profiles of steady-state c{x) (blue dots) and a{x) (dashed red) 
are shown in Figs. \V[h) and (c) from the smooth and discon- 
tinuous regions, using (r, A) = (0.9,10) and (1.2,15), re- 
spectively. For all figures presented in this Letter we use pa- 
rameters a = 7 = \/0A, /3 = 0, ^ = 1, and w = 0. For 
comparison, in Figs.fljd) and (e) we show c(n) and a{n) pro- 
files inferred from numerical simulation of the detailed ionic 
model in Ref. lITsll . where n indexes individual cells, using 
parameter values that give smooth traveling profiles and dis- 
continuous stationary profiles, respectively. Traveling profiles 
have arrows indicating movement. 

The onset of alternans at ri (A) is mediated by an abso- 
lute instability analogous to that studied in Ref. ^ for the 
voltage-driven case. For (3 = a linear stability analysis 
yields thresholds of ri(A) = I - rj + St]^'^^^ /{AA'^^^) and 
1 — ?7 + £,'^{wA)^^ for the instability of the traveling and sta- 
tionary modes, respectively, where rj = aj. Furthermore, the 
wavelength at onset is 47r^2/^A^/^/ V3 and 27r(ti;A)^/^ in the 
traveling and stationary cases, respectively, which agrees with 
the voltage-driven case in Ref. Similar expressions can be 
obtained for /3 7^ 0. Numerical simulations (not shown) are in 
good agreement with these theoretical results. 

We now concentrate on the discontinuous regime that is 
the primary focus of this letter To characterize calcium al- 
ternans profiles in this regime [cf. Fig. [lie)], we examine 
first stationary steady-state period-two profiles and substitute 
c(x) = Cn{x) = —Cn+i{x) into Eq. dHJ. After differentiating 
Eq. dH) with respect to x and some manipulations, we obtain 



Ac'(x) 



'(x) - (r - l)c{x) - aAa'{x) 
(r-l)-3c2(a;) ' 



(10) 



Thus, when alternans grow from c 



with r > r2(A), 



if c{x) — C- = ±^(r — l)/3 the derivative diverges and 
c{x) becomes discontinuous. Through the discontinuity, the 
quantity ad{x) in Eq. ^ remains smooth, so finding the 
other root of the cubic (r — l)c = c? + ad gives the value 
of c{x) at the latter end of the discontinuity. This gives 
c{x) = c+ = =f2-\/ (r — l)/3 and a total jump of ampHtude 
|c+ — c_ I = •\/3(r — 1). To measure the asymmetry at a node, 
we introduce the quantity A = ||c+| — |c_||/ \J {r — l)/3. 
We will refer to a discontinuity where c{x) jumps from c_ = 
(r — l)/3 to C-i- ~ =f2-\/ (?■ — l)/3 as a normal jump. A 
remarkable property of this jump is that the limiting values c+ 
and c_ on either side of the node depend only on the strength 
r of Ca; -driven instability, and is independent of all the other 
parameters A, 77, P, ^, and w. Experimentally, r = 1 is the 
point in parameter space where an isolated myocyte paced 
with a periodic AP-clamp waveform, or a tissue paced at one 
point with negligible CV-restitution (A = 00) bifurcates to al- 
ternans. Hence, the ratio c+ /c_ can be used to deduce r in 
tissue experiments or simulations under a finite effect of CV- 
restitution, and hence to relate single-cell and tissue behavior 
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FIG. 2: (Color online) Using r — 1.2, jc_| (solid blue) and |c+| 
(dashed red) versus A from an initial profile with A = 10. Inset: A. 

When starting from the unstable base solution without al- 
ternans (a = c = 0) in the regime r > r2{A), SDA forms 
dynamically as a periodic pattern of discontinuous nodes with 
normal jumps. A unique feature of SDA evolution in this 
regime, which is entirely absent for Kn-dominated instabil- 
ity, is that both the node positions and alternans profiles can 
depend strongly on the history of how the parameters r and 
A are varied. If A or r are increased starting from a profile 
with normal jumps, the position of the nodes remains con- 
stant, but the shape of the profile deforms in such a way that 
the jump in Cai-alternans profile becomes symmetrical about 
the node, i.e. both |c_| and |c+| approach the same limiting 
value \c±\ = \/r — 1 where A vanishes. This shows that, 
if initial conditions contain discontinuous nodes, jumps need 
not be normal in steady-state if c(x) = c_ = ±\/{r — l)/3 is 
not attained. This is shown in Fig. |2] where we plot |c_ | (solid 
blue) and |c+| (dashed red) using r = 1.2 as we increase A 
from 10. As A is increased, |c_ | and |c+ 1 tend toward one an- 
other We superimpose theoretical values of \c± \ for the nor- 
mal jump and A^^ = cases in dot-dashed black for compari- 
son, noting that | c_ | , | c+ 1 vary smoothly between these values 
for intermediate values of A. In the inset we plot A versus A, 
noting that A —> as A —> 00. If A is then decreased back un- 
til it reaches its original value (not shown), the profile recovers 
its original shape. However, if A or r are decreased starting 
at a point where the jumps are normal, the pattern close to 
the node preserves its shape, but the node moves towards the 
pacing site. Importantly, if A or r are increased back after the 
node has moved, the node does not return to its original posi- 
tion, but rather its shape will deform to become symmetrical 
as described above. Since no parameter change can induce 
the node to move away from the pacing site, node motion is 
unidirectionally pinned. We note that we have also observed 
unidirectional pinning in our ionic model simulations 11511 . 

When the node is unpinned, we find that the location of 
the first node, denoted xi, scales linearly with A, suggesting 
that the node spacing is independent of electrotonic coupling. 
This linear scaling with A in the discontinuous regime is to 
be contrasted with the scaling of the node spacing for smooth 
alternans profiles (e.g. ^2/3^^1/3 fQj- ^ _ q-j^ which depends 
strongly on electrotonic coupling. Physically, this linear scal- 
ing reflects the fact that electrotonic coupling has a negligi- 
ble effect on the outer scale where the alternans profile varies 
slowly on a scale ^ A, and only becomes relevant on a scale 
^ ^ near the nodes. This only adds a subdominant correction 
of order ^ to the xi ^ A scaling. Mathematically, it can be 
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FIG. 3: (Color online) (a) Paths (i) and (ii) (solid and dashed) in 
(r, A), (b) Initial profile (solid black) and final profiles (solid and 
dashed blue) after moving along path (i) and (ii). (c), (d) Asymmetry 
A (circles) and xi (crosses) along paths (i) and (ii), respectively. 

related to the fact that A scales out of Eq. ( fTOb in the limit 
A ^ if one uses the scaled variable x = x/ A instead of x. 

To investigate the consequences of unidirectional pinning, 
we investigated the pattern evolution in response to multi- 
ple parameter changes. Namely, we changed r and A fol- 
lowing two different paths that connect the same points in 
(r, A). Starting with a pattern with normal jumps at (r, A) = 
(1.16,30), we move to (1.26, 14) first by increasing r and 
then decreasing A [path (i)] and vice versa [path (ii)]. De- 
spite the same start and end parameters, the resulting profile 
characteristics vary significantly depending on the path fol- 
lowed, as shown in Fig. [3] In Fig. 13 a) paths (i) and (ii) are 
denoted by solid and dashed blue arrows, respectively, with 
the start and endpoints denoted as a black circle and square. 
In Fig. Ob) we zoom in on the first node of the initial pro- 
file (solid black curve) and final profiles after moving along 
path (i) (solid blue) and (ii) (dashed blue). Consistent with 
our summary above, Xi remains constant along path (i) while 
A decreases as r is increased, after which A increases as A is 
decreased. Along path (ii) xi decreases with A, then remains 
constant while A decreases as r is increased. This is shown in 
Fig- [3(c) and (d), which show A (black circles) and xi (red 
crosses) measured along paths (i) and (ii), respectively. 

In conclusion, we have extended our basic theoretical 
understanding of SDA dynamics to the important case of 
calcium-driven instability. Furthermore we have made a num- 
ber of new experimentally testable predictions, which we have 
validated by detailed ionic model simulations. The main pre- 
diction is that node motion becomes unidirectionally pinned 
when the Cai alternans profile becomes spatially discontinu- 
ous above a threshold of Caj-driven instability. This predic- 
tion could be tested by first increasing progressively the pac- 
ing rate (decreasing the inverse restitution slope A), thereby 
causing the node to move towards the pacing site, as has al- 
ready been observed in some experiments on APD SDA flS, 
and then decreasing the pacing rate to its original value. If 



the Ca^ -alternans profile exhibits a jump at the node, the 
node should remain stationary. Increasing pacing rate can 
also cause several parameters to change, including the de- 
gree of Caj-driven instability. However, we have shown that 
history-dependent SDA evolution is robust to multiple param- 
eter changes (Fig. O, and hence should be observable in more 
complex situations. We emphasize that unidirectional pin- 
ning is a purely dynamical phenomenon independent of in- 
trinsic tissue heterogeneities, which can also potentially pin 
node motion. However, we expect pinning due to tissue het- 
erogeneities to be generally bi-directional, and hence distin- 
guishable from unidirectional dynamically-induced pining. A 
second prediction is that the spatial jump in Cai-alternans am- 
plitude displays remarkably universal features. The magni- 
tude and asymmetry of this jump are insensitive to most pa- 
rameters except the degree of Cai-driven instability, and both 
quantities are generally history-dependent. 

Unidirectional pinning generally makes it harder to elimi- 
nate SDA by node motion once they are formed. We there- 
fore expect SDA to be more arrhythmogenic for Cai- than 
-dominated instability. Given that alternans are believed 
to be predominantly Caj-driven in common pathologies such 
as heart failure, SDA may play an even more important role 
than previously thought in such pathologies. 
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